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I. INTRODUCTION 

One arena where non-perturbative QCD could be 
tested in the future is in the study of nuclear matter 
in extreme conditions - high density or high tempera- 
ture. An approximation to this situation will eventually 
be produced in relativistic heavy ion collisions. In con- 
trast to the familiar confining phase of QCD which is 
relevant to nuclear physics and where the degrees of free- 
dom are hadrons, which are the colorless bound states 
of quarks and gluons, at sufficiently high temperature or 
density the system should be a plasma of the unconfined 
quarks and gluons themselves. Although there may not 
be a distinct phase transition between these regimes, the 
existence of the quark gluon plasma should have some 
well-defined signatures and studying its properties is im- 
portant. 

At zero temperature, for states near the vacuum of 
QCD, it is very clear what confinement means - there are 
no asymptotic states of the quantum field theory which 
carry color charge - all states are singlets of the global 
color symmetry. On the other hand, at finite temperature 
or density, where the physical system is in a mixture 
of states, there seems to be no clean characterization of 
confinement. Part of the problem is the fact that the 
confining interaction is so strong at large distances that it 
screens itself. If a single test quark carrying color charge 
is introduced into a confining medium, there is a low 
energy state where the color charge of the test quark 
is screened by a dynamical anti-quark to make a color 
neutral object with finite energy. 

Of course, in a hypothetical system where quarks are 
absent, or where they occur only in the adjoint represen- 
tation of the gauge group, so that they cannot screen the 
color charge of a fundamental representation test quark, 
one can characterize confinement by asking how much en- 
ergy it takes to introduce the test quark into the system. 
In the confining phase it would be expected that it takes 
infinite energy, in the deconfined phase, finite energy. 



In the following, we will review a study of the issue of 
confinement in a particular high temperature gauge the- 
ory which is a toy model of QCD. The idea is to study 
a model which is simple enough to be exactly solvable, 
but still complex enough to exhibit the phenomenon of 
interest - a phase transition between quark confining and 
a plasma phases. The result of this study will be some 
speculations about characterizing the confining phase of a 
gauge theory at finite temperature and density when fun- 
damental representation quarks are present. The model 
with adjoint representation quarks can be solved in any 
dimensions 0. However, here, we will review only the 
one-dimensional case M . 

In two dimensions, the Yang-Mills field itself has no 
propagating degrees of freedom. In adjoint QCD, the 
matter fields provide dynamics by playing a role analo- 
gous to the transverse gluons of higher dimensional gauge 
theory. In fact, dimensional reduction of three dimen- 
sional Yang Mills theory produces two dimensional QCD 
with massless adjoint scalar quarks. Moreover, since ad- 
joint matter fields do not decouple in the infinite N limit, 
the large N expansion is of a similar level of complexity 
to that of higher dimensional Yang-Mills theory. One 
would expect it to exhibit some of the stringy features of 
the confining phase which are emphasized in that limit. 

Although adjoint QCD is not explicitly solvable, even 
at infinite N, details of its spectrum were readily an- 
alyzed by approximate and numerical techniques |5]-[|. 
In addition, Kutasov Q exploited an argument which 
was originally due to Polchinski |t]] to show that the con- 
fining phase must be unstable at high temperature and 
suggested it as a tractable model where the confinement- 
deconfinement transition could be investigated. 

Here, we discuss a simplified version of 1+1 dimen- 
sional QCD.f] This model was formulated in ref. Q and 



'Generalization of these ideas to higher dimensions is pos- 
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jjj . We consider a one-dimensional gas of non-dynamical 
particles which have color charges and which interact 
with each other through non-Abelian electric fields. 

Even when the quarks are in the adjoint representa- 
tion, this system exhibits confinement. Because there are 
no dynamical gluons which could screen adjoint charges 
in one dimension, at low temperature and density, ad- 
joint quarks are confined]^ They form colorless "hadron" 
bound states with two or more adjoint quarks connected 
by non-dynamical strings of electric flux. The large N 
limit resembles a non-interacting string theory in that 
the energy of a state is proportional to the total length of 
all strings of electric flux plus a chemical potential times 
the total number of quarks. The property of confinement 
is defined by estimating the energy required to introduce 
an external fundamental representation quark-antiquark 
pair into the system. In the confining phase, where the 
hadron gas is dilute, the quark-antiquark energy is pro- 
portional to the length of the electric flux string which, to 
obtain gauge invariance, must connect them. This gives 
them a confining interaction which increases linearly with 
their separation. 

In the confined phase, the average particle number den- 
sity and the energy density are small — in the large N 
limit both are of order one, rather than iV 2 which one 
would expect from naive counting of the degrees of free- 
dom. This is consistent with the fact that in a confining 
phase the number of degrees of freedom, i.e. hadrons, is 
independent of N. In contrast, in the deconfined phase, 
since the number of degrees of freedom, i.e. quarks and 
gluons, is proportional to iV 2 the particle density and 
energy are also of order N 2 . 

As temperature and density are increased, eventu- 
ally we arrive at the situation where inserting a quark- 
antiquark pair, involves a negligibly small addition to the 
energy of the total system. This is the signature of the 
deconfined phase. 

Between these two phases is a transition, which we 
show, is of first order B. In the string picture, this phase 
transition occurs when the strings in a typical configura- 
tion percolate in the one-dimensional space. The order 
parameter is the Polyakov loop operator ||,[l0| which 
measures the exponential of the negative of the free en- 
ergy which is required to insert a single, unpaired fun- 
damental representation quark source into the system. 
This free energy is infinite (and the expectation value of 
the Polyakov loop is zero) in the confining phase and it 
is finite in the deconfined phase. 



sible and a number of models are explicitly solvable there as 
well El. Here, for simplicity we concentrate on the one di- 
mensional case. 

^In higher dimensions, an adjoint charge and a gluon could 
form a color singlet bound state. 



The formalism we will use allows for a straightforward 
extension of the model to include non-dynamical, funda- 
mental representation 'quarks' j^j in addition to adjoint 
representation 'gluons'. QCD with fundamental repre- 
sentation quarks is solvable in the large N limit and 
has been used as an explicit example where confinement- 
related phenomena can be studied |llj. Its solvability 
follows partially from the fact that the quarks decou- 
ple in the infinite N limit. We can also solve the non- 
Abelian coulomb gas when it has fundamental represen- 
tation quarks, which couple in a non-trivial way when 
their fugacity is of order N. This resembles a field the- 
ory when the number of flavors of heavy quarks, Np, is 
of the same order as the number of colors. Then our 
solution applies to the large N, Np ~ N, heavy quark 
limit. 

In a gas dominated by fundamental representation par- 
ticles, the energy needed to introduce another fundamen- 
tal representation quark is always finite, and more so- 
phisticated ideas are needed to study confinement p2| . 
(For an interesting suggestion about symmetry breaking 
in this case, see ref. Intuitively this is because the 

electric flux associated with any external fundamental 
representation source can be screened by a fundamen- 
tal representation quark which is already available in the 
system. We will find that there is a third order phase 
transition in this system, where the character of the con- 
fining phase changes in a fundamental way. 

The system we will consider here is that of a mixed gas 
of adjoint and fundamental representation non-Abelian 
charges interacting via 1+1 dimensional U(N — > oo) 
gauge fields. The main result of our investigation is the 
development of a unified order parameter [ fl4| for discern- 
ing the low temperature confined phase from the high 
temperature quark-gluon plasma phase for all relative 
densities of quarks and gluons. This order parameter can 
be defined by considering the group theoretic details of 
the mean-field solution of the model in the large N limit. 
We will outline this formalism and give some examples 
of its utility. 

II. FORMALISM 

The action of 1+1-dimensional Yang-Mills theory cou- 
pled to a number of static particles is 




In the canonical quantization of this system, the electric 
field is the canonical conjugate of the spatial component 
of the gauge field, 
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[A a (x),E b (y)] = i5 ah 5{x-y) 
The Hamiltonian is^] 

„2 N 2 



J 0=1 



(2) 



(3) 



and the temporal component of the gauge field, Aq plays 
the role of a Lagrange multiplier field which enforces the 
Gauss' law constraint, 



^-E a (x) ~ f abc A b (x)E c (x) 
ax 

K 



(4) 



i=l 



We will impose this constraint as a physical state con- 
dition. The particles with color charges are located at 
positions Xi, . . . , xk- are generators in the represen- 
tation Ri operating on the color degrees of freedom of 
the i'th particle. 

In the functional Schrodinger picture, the states are 
functionals of the gauge field and the electric field is the 
functional derivative operator 



E a (x) 



1 



i 8A a (x) 



The functional Schrodinger equation is that of a free par- 
ticle 



dx — — 



A;xi,.. .,x K ] 
S 2 



e 2 ^ ^ 



2 ^ (SA-(xW 



(5) 

[A;xi, . . .,x K ] 



Gauss' law implies that the physical states, i.e. those 
which obey the gauge constraint (Q), transform as 

* ai " a «[A B ;x u ...,x K ] (6) 
= 9aX • ■ ■ & K {x K )^- hK [A; x 1; ... , x K ] 

where A 9 = gAg^ + igVg^ is the gauge transform of A. 
For a fixed number of external charges, this model is ex- 
plicitly solvable. In the following we shall examine its 
thermodynamic features, where we assume that the par- 
ticles have Maxwell-Boltzmann statistics. 

We shall work with the grand canonical ensemble. The 
partition function is constructed by taking the trace of 
the Gibbs density e~ H l T over physical states. This can 



"'•Here, for concreteness, we consider U(N) gauge theory. The 
gauge field is A = A a t a , with t a the generators in the funda- 
mental representation. 



be implemented by considering eigenstates of A a (x) (and 
an appropriate basis for the non-dynamical particles) 



\ A )e ai 



Projection onto gauge invariant states in- 



volves a projection operator which has the net effect of 
gauge transforming the state field at one side of the trace, 
and integrating over all gauge transformations |D| . The 
resulting partition function is 

Z[xi,T]= (7) 
= J[dA][dg] (A\e- H ' T \A 9 )Tr g^{ Xl ) . . . Tr 9 Rk (x k ) 

where [dg(x)] is the Haar measure on the space of map- 
pings from the line to the group manifold and [dA] is 
a measure on the convex Euclidean space of gauge field 
configurations. Here, we will consider the case of par- 
ticles in both the adjoint and fundamental representa- 
tions. For U(N), the trace in the adjoint representation 
is Tr g Ad (x) — |Tr g(x)\ where g(x) is in the fundamen- 
tal representation. In order to form the grand canoni- 
cal ensemble, we average over the particle positions by 
integrating J dx\ . . . J dxK, multiply by the fugacity, A 
for adjoint charges, K for fundamental charges, k* for 
the conjugate to the fundamental representation, to the 
power of the number of respective charges, divide by the 
factorial statistics factor to obtain Maxwell-Boltzmann 
statistics and sum over all numbers of particles. The 
result is 



Z[X,k,T] 



[dA][dg] 



-S etf [A,g] 



(8) 



where the effective action is 



-S B{i [A,g] 



(A 



-H/T 



A 9 ) 



(9) 



x exp (^J dx (A |Tr g\ 2 + kTt g + K*Tr 



The Hamiltonian is the Laplacian on the space of gauge 
fields. Using the explicit form of the heat kernel 



{A 



-H/T 



A 9 ) ~ exp - 



[ dx J Tr (A - A 9 ) 



we see that the effective theory is the gauged principal 
chiral model with a quadratic potential 

S cS [A,g}= (10) 

J dx (J^- Tr \Vg + i[A, g}\ 2 - A |Tr g\ 2 - 2Nk ReTr g^j 

Here the integration over gauge fields A effectively en- 
forces Gauss' law as one integrates over all elements of 
the gauge group with the Haar measure [dg] . The fugac- 
ities of the adjoint and fundamental charges are given by 
the parameters A and Nk, respectively. Since we con- 
sider the matrix-valued fields A and g to be taken in the 
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fundamental representation of SU(N), the large N limit 
will lead directly to the familiar situation of matrix mod- 
els with large N x N matrices. In order to keep all terms 
in the action of ( |ic| ) at leading, TV 2 , order in this limit we 
will restrict parameters of the system such that 7 = -?m , 
A and k are each of 0(1). 

For the discussion of a confinement-deconfinement 
phase transition the most important aspect of the action 
in ([h]) is a global symmetry S[A, g] = S[A, z g] when the 
fundamental charge fugacity n vanishes. Here z is a con- 
stant element from the center of the gauge group, which 
for U(N) is U(l) and for SU(N) is Z N . It is this symme- 
try and its (thermo-)dynamical breaking that leads to the 
deconfinement phase transition in this model. If k ^ 
the question of what remnants of this symmetry persist 
is one we will answer in the next sections. 

Additionally, there is a gauge invariance that can be 
used to diagonalize the matrices gij(x) — e loli ^5ij. The 
density of eigenvalues p{9,x) = J2iLi ~ a i( x )) 
corresponding to the large N saddle-point evaluation of 
( |l0| ) now completely characterizes the properties of the 
system. Our goal is to find this distribution of eigen- 
values. Without loss of generality we can consider the 
Fourier expansion 



Zir Zn * — ' 



rie ~ me , c n (x)* = c- n (x) 
(11) 



The configurations of the eigenvalue density (|ll|) that 
saturate © at large N can be found via the collective 
field theory approach |l6|JT7| ]. The method is essentially 
based on the relation between matrix quantum mechanics 
and non-relativistic fermions Jll|[nJ . Leaving the details 
to H , it can be shown that a solution of the saddle-point 
evaluation of fllO| ) is given by 



Po(d) 



■^\JE + 2(Aci -(- k) cosB where p is real 



otherwise 



(12) 



The constant of integration E has physical interpretation 
as the Fermi energy of a collection of N fermions p8[ on 
the circle subject to a periodic potential. It is fixed by 
requiring the eigenvalue distribution to have unit normal- 
ization. Furthermore, since the potential due to adjoint 
charges is non-local in eigenvalue space, the Fourier mo- 
ment c\ (see (pT|)) must be self-consistently determined 



Cl 



d6 pq{9) cos( 



(13) 



This pair of conditions is most conveniently analyzed by 
introducing a new parameter p = E/(2(\c± + and the 
integrals over the positive support of p + cos 8, 




Kly 

FIG. 1. Plot of the lines (2.7) for p ranging from 0.4 
(upper right corner) to 75 (line at the extreme left). The 
region of overlapping lines corresponds to a region of first 
order phase transition. 



In{p) = — [ d6 cos nO^J p + cos ( 

7T J 



(14) 



In terms of p, the solution of the normalization and mo- 
ment conditions is given by 



('1 



h(p) 



and 



K 

1 



1 



A IiQx) 

4/ (/i) 2 7 / (/i) 



(15) 



(16) 



This last relation gives a family of lines in the (A/7, ft/7) 
plane parameterized by p. As is shown in Figure this 
family overlaps itself for lower densities of fundamental 
charges, /t/ 7 signaling the fact that there are multiple 
solutions to the equations of motion in this region of the 
phase diagram. Considering the free energy one can show 
Q that for lower densities of adjoint charges the stable 
solution has p > 1 while at higher adjoint densities the 
stable solution has p < 1. In the intermediate regime 
there is a first order phase transition. As the density of 
fundamental charges is increased the first order transi- 
tion is smoothed out and a third order phase transition 
persists along the line p = 1. 

The parameter p is now seen to be useful for two differ- 
ent reasons. First it characterizes the general structure of 
the phase diagram (Figure ^|) where the 'strong coupling' 
regime is the region with p > 1 and the 'weak coupling' 
regime has p < 1. As well, and of more importance for 
our analysis, we find that the expectation values of traces 
of powers of the group element g are given as a function 
of the single parameter p 



irn n / Art In{p) 
( ^ 9 /N) = n = ~h(p) 



(17) 



Consequently, it makes sense for our purposes to re-define 
the eigenvalue distribution in terms of p 
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FIG. 2. Schematic picture of the phase diagram for the ad- 
joint and fundamental representation non-Abelian Coulomb 
gas. The dotted curve marks the first order part of the critical 
line. The solid curves above and below it are the boundaries 
of the area with two possible phases. They join at a point 
which shows second order behaviour. For larger /t/7, we find 
a third order line (fi = 1) marked by a solid line. 



po(o, m) = —rr\ Vv + cos0 (18) 

ttIq (fx) 

In the next section we will use this definition and its 
connection to the dominant configuration of the gauge 
element to analyze the phase diagram in terms of group 
theory. 

III. ORDER PARAMETERS AND THE THEORY 
OF GROUP CHARACTERS 

As is known, in the case of pure gluo-dynamics, the 
realization of the center symmetry of the gauge group 
governs confinement p]lCf|. The Polyakov loop op- 
erator Tr g(x), which is related to the free energy 
— Tlog (Trg(x)Trg 1 '(0)) of a conjugate pair of static, ex- 
ternal fundamental charges separated by a distance x, 
serves as an order parameter f2lj| to test confinement. 
Since Tr g(x) transforms under the center as Tr g(x) — > 
z Tr g(x), the expectation value of the Polyakov loop 
operator must average to zero if the center symmetry 
is preserved. Physically this suggests that an infinite 
amount of energy is required to introduce a single fun- 
damental test charge into the system. The presence of a 
gas of fundamental charges (n ^ 0) changes this situa- 
tion though by explicitly breaking the center symmetry. 
Consequently we lose the Polyakov loop operator as an 
order parameter for phase transitions in the system. In 
this section we introduce a suitable generalization of the 



Polyakov loop operator which will allow us to identify a 
new order parameter. 

As seen in the previous section, the solution of the 
non-Abelian Coulomb gas with adjoint and fundamental 
representation charges is completely characterized by a 
Fourier sum of the traces c n = (Tr g n /N) - the higher 
winding Polyakov loops. As noted in j| the character 
of these traces changes between the strong and weak 
coupling regimes. In particular, in the strong coupling 
(// > 1) phase c„ is damped exponentially with n while 
in the weak coupling (/i < 1) phase the damping follows 
a power law behaviour. Now we will reconsider this be- 
haviour in terms of group theory. 

Since the matrix g is an element of the special unitary 
group, its trace in an irreducible representation, R defines 
the group character for that representation 

Xfl(s) = Tr fl g (1) 

For the N dimensional fundamental representation of 
SU(N), F, the group character is just the Polyakov loop 
operator described above since we are considering group 
elements to be taken in the lowest fundamental represen- 
tation 

XF (g) = Trg (2) 

Further simple examples are the symmetric (S) and anti- 
symmetric (A) combinations of a pair of fundamentals 
where we have 

Xs ( 5 ) = l[(Tr 5 ) 2 + Tr ff 2 ] (3) 

and 

^(. 9 ) = i[(Tr. 9 ) 2 -Tr 5 2 ] (4) 

A general relation between characters and the group ele- 
ments is given by the Weyl formula but is not necessary 
for the following. A complete discussion can be found in 
standard references (see |22| for example). 

The main idea is that the eigenvalues of the group ma- 
trices, which are the only relevant dynamical variables 
in the grand partition function (|l0|), are completely de- 
termined by the N quantities {Tr g n }, n = 1 ... N. In 
turn these traces form an algebraic basis equivalent to 
the characters of the N fundamental (completely anti- 
symmetric) irreducible representations of SU(N) (in- 
cluding the trivial representation) . Here we will explicitly 
demonstrate the relationship between the basis of traces 
and the basis of group characters. Ultimately it is the 
group theoretic variables which we will use to character- 
ize the phases of the model ([to]) . 

The standard basis for general functions (of finite de- 
gree) of the eigenvalues of a matrix is the set of elemen- 
tary symmetric functions {a r }. In terms of the eigen- 
values Xj — e %6 i of the group element g they are given 
by 
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j 

j<k 

as = 22 ^j^kM 

j<k<l 



a at = J^J A j = dot .9 = 1 



(5) 



(G) 



with a,. = for r > N. The relationship of the symmetric 
functions {a r } to the traces of the group elements, S n = 
Tr g n , is given p3|] by the determinant 



ak= k\ 



51 1 

5 2 5i 2 

53 S*2 5i 

Sk-2 Sfc-3 

Sk Sk-i Sk-2 





£2 Si k — 1 
S3 5*2 Si 



(7) 



Most importantly, it can be shown that the elementary 
symmetric functions are nothing more than the charac- 
ters of the fundamental representations for the unitary 
group ^IH). That is, for the fundamental representa- 
tion which is the anti-symmetric combination of fc, N 
dimensional representations, Xk{d) — a k- 

The determinant ^ can be evaluated [pfl in terms of 
a multinomial expansion most compactly stated in terms 
of a generating function 



(8) 



z=0 



For our purposes though it is useful to convert to a con- 
tour integral about the origin. 



Xk(g) = 



(-1)* 

2ni 



dz 



y> Trg r ' 



I) 



(9) 



These last two expressions explicitly demonstrate the re- 
lationship between the group element g and the k th fun- 
damental representation of the gauge group and are com- 
pletely general results. 

With these relations we see that there is a direct con- 
nection between the gauge group element g and the 
irreducible (fundamental) representations of the gauge 
group. In particular, in the previous section we have 
seen that in the large N solution of the non-Abelian 
Coulomb gas a certain configuration of the gauge ma- 
trix, go saturated the evaluation of the partition function 
(|lTJ). Now it is natural to ask what is the configuration 
of irreducible representations corresponding to the dom- 
inant go . This corresponds to evaluating the expecta- 
tion (xk{g)) in the background of the non-Abelian gas. 



In principle this involves calculating expectations of the 
form (Tr g ni ■ ■ • Tr g rLr ) but because of the factorization 
of gauge invariant objects in the limit N 00, this re- 
duces to a product of expectations, (Tr g ni ) ■ ■ ■ (Tr g nr ). 
Consequently (Xfe(<?)) is determined by replacing Tr g n 
by its expectation value in (|^). Of course expectation 
values of the group element traces are intimately related 
to the eigenvalue density p(9,n) (see ([TT|)) hence, after 
performing an infinite sum, we obtain 



dz 

— exp 



(-!)« 



N 



2m 

d9p(6, rflog 



(10) 



2z cos 9 



,2a 



Note that we have defined a new real parameter a = k /N 
on the unit interval that effectively labels the fundamen- 
tal representations in the large N limit. Of course ( |l0| ) 
now depends on a continuous variable and is of a slightly 
different functional form than the discrete case (xfe)- I n 
the remainder of this discussion we will consider only the 
character parameterized by a as defined in (JToj) . 



IV. CALCULATION OF THE EXPECTATION OF 
FUNDAMENTALS 

In this section we will concentrate on calculating 
(Xa) with eigenvalue density (|l8|) for the non-Abelian 
Coulomb gas. This calculation will give a clear picture 
of the group theoretic excitations present in different re- 
gions of the phase diagram and consequently allow us to 
define an order parameter for the deconfinement transi- 
tion even in the presence of fundamental matter. 

Since explicit evaluation of (|l0|) is difficult we begin 
with some special limiting cases. As p — > — 1 the support 
of the eigenvalue distribution ( fL8[ ) vanishes at 6 = 0. 
The distribution does not vanish though as it retains unit 
normalization and effectively becomes a delta function, 
S(9). Consequently we find the gauge matrix g is just the 
identity at p = — 1, hence 



(Xa) = Jim 



N 
aN 



= 2 



N 



Nir 



-27V(a-l/2) 2 



(1) 



In this limit we find that the distribution of characters 
is symmetric about a = 1/2 as one would expect in a 
system where the total colour charge is vanishing. As well 
in this limit (xa) is non- vanishing and all fundamental 
representations are present in the large N background 
solution of the model. As we will see, this result is generic 
in the weak coupling phase p < 1. 

In the opposite limit, as p — > 00, it can be shown 
that the eigenvalue distribution ( [l8| ) approaches a con- 
stant value p = X/2tt with the eigenvalues of the group 
element g becoming uniformly distributed on the unit 
circle. Since expectation values of the traces of powers of 



the gauge matrix are essentially Fourier transforms of the 
eigenvalue distribution, it is easy to see that (Tr g n ) — > 
in this limit and 



(Xa) -* $0,a 



(2) 



This limit corresponds to the extreme strong coupling 
phase of the model where the Polyakov loop operator 
((Tr g) r~j (xi/n)) has vanishing expectation value and 
the standard analysis would point to a phase where colour 
charges are strictly confined into hadron-like structures. 

In general the integral (|i"o| ) can be evaluated by saddle- 
point methods in the large N limit in which we are in- 
terested. The relevant action in this limit is 



S(a,n, z) = / d9p(9, /i) log 



1 



2z cost 



(3) 



Solving the stationarity condition, dS/dz\ Zo = 0, for a 
in terms of zq we find the saddle-point condition for the 
large N behaviour of the integral (|l^) is given by the 
relationship 



d6p(6, n) 



z o( z o 



cost 



1 + Zq — 2z cos 9 



(4) 



Since a is a real parameter restricted to the unit interval 
[0, 1] it can be shown that the saddle-point value of the 
parameter zq is real. Further, for zq > 1 and < zq < 1 
EqnJJ returns values of a > 1 and a < 0, respectively. 
Consequently we need only consider real, negative values 
of the parameter Zq. 

We now turn to an examination of the saddle-point 
approximation of ( |i"o| ) for different regions of the phase 
diagram of the model at hand beginning with the weak 
coupling phase, /i < 1. In this case the support of 
the eigenvalue distribution ( |i~8| ) is bounded away from 
9 = ±7r and hence the denominator in (Q) is non-singular 
for all values of Zq. Consequently, in this regime a varies 
smoothly and monotonically with zq and the relation (Q) 
can in principle be inverted to obtain zq(o). With this 
information, the large N asymptotic form of the expec- 
tation value of the characters (xa) can be determined by 
standard saddle-point methods. In Figure [| we show a 
numerically calculated example of a as a function of zq 
for fi — 0.5. For this same case we show a schematic dia- 
gram of the magnitude of the expectation value | (Xa) I as 
a function of a in Figure ^. In particular we see that the 
system has excitations in all irreducible representations. 

For /i > 1 the situation is somewhat different. Now 
the support of the eigenvalue distribution is the full 
interval 9 S [— n, it], and the denominator of (Q) causes 
non-analytic behaviour to appear. As one increases fi 
through unity the saddle-point relation for a shows this 
non-analytic behaviour as a discontinuity at zq = — 1 (see 
Figure |5|). The result is that an open interval of a values 
centered on a = 1/2 are mapped into this discontinu- 
ity when the saddle-point relation (H) is inverted. Since 
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FIG. 3. Plot of the solutions of the saddle-point relation 
(7.60) for n = 0.5. 



I<x a >l 




FIG. 4. Schematic diagram of |(xa)| vs - a f° r M = 0-5- 
Note that all fundamental representations have non- vanishing 
expectation value. 
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becomes too energetically costly to produce and the sys- 
tem can no longer screen the interaction between a pair 
of a — 1/2 fundamental charges. In this strong coupling 
phase the interacting pair sees a linear confining poten- 
tial (though somewhat reduced as compared to the empty 
background). As one further increases fi the gap in the 
spectrum of fundamental representations becomes larger 
and in the extreme limit /i = oo the system contains only 
excitations in the trivial representation. This is precisely 
the confining phase of pure gluo-dynamics. 
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FIG. 5. Plot of the solutions of the saddle-point relation 
(7.60) for n = 1.2. 




0.6 

a 



FIG. 6. Schematic diagram of |(x<*)l vs - a f° r M = 1-2- 
In this case the expectation value of representations with a 
between ~ 0.25 and ~ 0.75 is vanishing. 



this discontinuity occurs in the saddle-point relation, it is 
not surprising to find that the curvature associated with 
the Gaussian integration of the saddle-point approxima- 
tion is divergent, effectively forcing the integral to vanish. 
In terms of the expectation values of different represen- 
tations in the background of the non-Abelian Coulomb 
gas, we see that an open interval of fundamental repre- 
sentations centered about a = 1/2 is missing from the 
spectrum in the large N limit. In Figure ^| we show an 
example of the behaviour of the expectation value |(Xa)| 
with a for fi = 1.2. 

The main outcome of this analysis is that the expec- 
tation value of the central fundamental character (xi/2) 
is vanishing if and only if fi > 1. Consequently it may 
be considered an order parameter distinguishing between 
the strong and weak coupling phases of the model. Phys- 
ically the situation is clear. In the weak coupling phase 
the system can effectively screen the interactions of any 
pair of charges regardless of their representation since 
the system contains excitations in all representations of 
the gauge group. We conclude that the system looks 
much like a quark-gluon plasma where charges are ef- 
fectively deconfined. At the phase transition line non- 
Abelian flux in the a = 1/2 fundamental representation 
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